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ABSTRACT 

We show, through a simple patchy reconnection model, that retracting reconnected 
flux tubes may present elongated regions relatively devoid of plasma, as well as long 
lasting, dense central hot regions. Reconnection is assumed to happen in a small patch 
across a Syrovatskii (non-uniform) current sheet (CS) with skewed magnetic fields. The 
background magnetic pressure has its maximum at the center of the CS plane, and 
decreases toward the edges of the plane. The reconnection patch creates two V-shaped 
reconnected tubes that shorten as they retract in opposite directions, due to magnetic 
tension. One of them moves upward toward the top edge of the CS, and the other one 
moves downward toward the top of the underlying arcade. Rotational discontinuities 
(RDs) propagate along the legs of the tubes and generate parallel super-sonic flows that 
collide at the center of the tube. There, gas dynamics shocks that compress and heat the 
plasma are launched outwardly. The descending tube moves through the bottom part 
of the CS where it expands laterally in response to the background magnetic pressure. 
This effect may decrease plasma density by 30 % to 50 % of background levels. This 
tube will arrive at the top of the arcade that will slow it down to a stop. Here, the 
perpendicular dynamics is halted, but the parallel dynamics continues along its legs; 
the RDs are shut down, and the gas is rarified to even lower densities. The hot post- 
shock regions continue evolving, determining a long lasting hot region on top of the 
arcade. We provide an observational method based on total emission measure and 
mean temperature, that indicates where in the CS the tube has been reconnected. 

Subject headings: magnetic fields — magnetohydrodynamics (MHD) — shock waves — 
Sun: flares 



Introduction 



Reconnection, the mechanism where field lines are broken and connected to other ones with 
different topology, is present in many theories that undertake to explain the initiati on and evolu- 
tion of the biggest explosions in the solar syst e m, flares or co ronal mass ejections (ILin fc Forbes 



2000 



Forbes fc ActonTll996l : Forbes et alJbood : ISchriiverlbood ): see also (jForbesl l200oh for 



energy 
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estimations in coronal mass ejections (CMEs)). Large amounts of stored magnetic energy can be 
released as thermal, ki netic, and radiat ive energy, as well as accelerated particles, when field lines 
change their topology (jMandrinil I201Q ). Reconnection is also involved in the interaction of the 
solar wind with the Earth's magnetosphere, where substorms that may affect our space-weather, 
electronics and telecommunications, are originated. For a recent review on magnetic reconnection 
in space plasmas and laboratory, see lYamada et al.1 (|201(J ) . 



Direct observation of the actual reconnection process has been e lusive. Its intern al structure 
remains unknown, even for well-diagnosed laboratory experiments (jKatz et al.l |2010| ). In order 
to explain the short times involved in the mentioned astrophysical explosions, fast reconnection 
(jPetschekl Il964l ) ought to occur in a minuscule diffusion region on a CS where non-ideal field-line 
transport coefficients are locally enhanced. This requires prohibitively high resolution measure- 
ments. Several fast reconnection models include localized magnetic resistivity enhancement on 
the CS rtllgai Tsudalll977l : IScholer fc RotiJ fl987l : ferkaev et alJbood : Ika fe Bhattacharieell200ll : 



Biskamp &: Schwarz 



2001 



All standard theoretical reconnection scenarios predict Alfvenic outflow s originated at the 
reconnection region. In steady state models, these outflows form coherent jets (|Parkerlll957l ; ISweet 



1958 



Petschek 



1964) , and in time-dependent models they are described as a localized retracti on of 



field lines (jSemenov et al.lll983l : iBiernat et alJll987l ; iHeyn fc Semenovlll996l : iNitta et al 



2001 



)• 



Only in the last fifteen years or so, has technology started to close the gap between theory and 
observations. Nonetheless, only indirect measurements of sola r corona reconnection have been pos 
sible . For example, the cusp-shaped, hot and dense flare loops (ITsuneta et al.lll 992; 



1996) consistent with the classical CSHKP model ( Carmichael 



(196 41): ISturrockl (|1968l ): Iffiravama 



Forbes h Acton 



(|l974l ) ; iKopp &: Pneumanl (jl976l ) ; see Figures 1 and 2 in iForbes &: Actonl (jl996l ) for a modern ver- 
sion including thermal conduction) are considered observational evidence of reconnection. These 
arcades, slowly growing and whose footpoints separate as their height rises, are interpreted as closed 
magnetic loops formed after two-dimensional, steady state reconnection. They pile up on top of 
each other to form the arcades. The earlier loops cool down as new loops lay on top of them, and 
an apparent footpoint motion is observed. The loops are hot, as expected from the conversion of 
magnetic energy into thermal energy from reconnection. 



There are several other indi r ect obseryation s of reconnection, like reconnection inflows (|Yokoyama et al 



2001 



Narukage &; Shibatall2006l : iLin et al.ll2005l ) with speeds between a few km s~* to approximately 



106 k m s" 1 , as well as Alfvenic (460— 3500 km s" 1 ) reconnectio n outflows ( Lin et al.ll2005l : IWang et al 



20071 ) and rec onnection in /out 



reconnection. 



Forbes Sz Acton 



1991 



) aboard Yohkoh (jOgawara et al. 




Sheeley fc Wangjl2007l ). Loop shrinkage also se ems to indicate 



using images of flare loops taken with SXT (ITsuneta et al 



1991 



), showed reconnected cusped loops shrinking to form 
post-flare arcades. This process takes hours and is different from the fast passage (1-2 minutes) 
of field lines through the reconnection outflows (jets). Hard X-ray (HXR) sources near the apex 
of soft X-ray loops indicate some high-energy process, possibly electron acceleration in that re- 
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gion (jMasuda et al.l 11994 ; iKrucker fc Linl 12008 ; iTomczakl |2009| ) . The acceleration may occur at 
the reconnection site or can be the result of the interaction of the downward reconnection outflow 



jet with the top of the SXR arcade, producing sup erhot plasmas (jShibata et al.lll995l ). Plasmoid 
(magnetic island) for mation (IForbes Priestl I1983T). as a consequ ence of reconnection, has been 
widely observed (e.g. IShibata et al.l (|1995l ); iNishizuka et al.l (I2010T)'). For a comp arative study of 
plasmoid in reconnecting CSs (solar and terrestrial contexts), see iLin et al.l (2003). 



Even though most of the mentioned observations can be understood by standard, steady state 
and two-dimensional reconnection theories, there are newly observed phenomena that cannot be 
reconciled with these representations. Supra-arcade downflows (SADs) and supra-arcade down- 
flowing loops (SADLs) seem to be two different observational s ignatures of retrac ting, isolated 
reconnected flux tubes with purely three-dimensional geometries (jSavage et al.ll2010l ). 



SADs were first reported more than one decade ago (jMcKenzie & Hudson 199 9]), and have 



been continuously observ ed since then by several instruments ( McKenziel 



200ll : llnnes et al.l l2003bl la): ksai et al.l hooi iMcKenzie Savage 



2009 



2000 



McKenzie Sz Hudson 



Savage et al.l l2010h . They 



descend through the supra-arcade fan (see ISvestka et al.l (|1998l ) for a descriptio n of one of such 
fans) toward the underly ing flare arcade and decelerate close to the arcade apex (jSheelev k, Wang 



2002 



Sheeley et al.ll2004l ). They are spatially localized as dark p ockets, followed by a dark wake or 



lane. Figure U shows ex amples of SADs, as observed by TRACE (IHandy et al 



1999 



2002 X-class flare (IWang et al 



1999 



Schrijver et al 



Golub et al.lll999) in 195 A pas s band. The image to the left shows an arca de for the April 21 



2002 



Gallagher et al 



2002 



Innes et al.ll2003al ). where two SADs 



are enclosed by a white dashed rectangle. Usually, arcades like this one can be described as two- 
dimensional, suggesting the presence of a wide spread CS above them. The images to the right 
are time sequences for the location delimited by the dashed rectangle. The white arrow on the top 
right snapshot indicates a dark pocket (SAD's head) that moves toward the arcade. This pocket is 
followed by a dark wake. The white arrow at time 1:44:33 indicates a second dark pocket (SAD) 
moving in the same direction. 

Downflows are observed as high as 40 — 60 Mm above the top o f soft X-ray arcades and 
interpreted as magnetic reconnection localized outflows (|McKenzid l2000l ) . Therefore, SADs appear 
to be the result of intermittent reconnection high in the corona. Three-dimensional and time- 
dependent models are needed to explain them. 

A three-dimensional model of patchy reconnection, accompanie d by magnetohydro dynamic 
(MHD) simulations, has been presented by Linton Sz Longcope (j2006l ). In that study, reconnection 
occurs across a one-dimensional, uniform CS with skewed magnetic fields. Resistivity is briefly 
enhanced by an unspecified mechanism within a sphere (reconnection patch) . This creates a pair of 
bent, thin reconnected flux tubes that retract in opposite directions. RDs prop agate along the legs 
of the tubes. Only the perpendicular dynamics of the tubes was analyzed by iLinton fc Longcope 
(|2006l ). who also argued that the descending coronal voids are reconnected flux tubes descending 
from a flare site high in the corona. 
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Longcope et abj (|2009i ) extended that study to include tube's parallel dynamics. The initial 
sharp angle at the reconnection site creates bends that move along the legs of the tubes as RDs. 
These bends rotate the magnetic field and generate parallel super-sonic inflows toward the center of 
the tube. There, the collision of these inflows launches outwardly two gas-dynamics shocks (GDSs) 
that heat and compress the plasma. The GDSs are disconnected from the resistive region, and are 
features of the ideal relaxation following reconnection. The rate of reconnection is not relevant in 
this case since the diffusion region is short-lived and does not participate in the post-reconnection 
evolution. 



In a recent paper (IGuidoni &; Longcopdl2010l ). we presented generalized thin flux tube (TFT) 
equations for the dynamics of reconnected flux tubes. These include a pressure gradient in the 
parallel direction, as well as temperature-dependent, anisotropic viscosity and thermal conductivity 
in non-ideal terms. In that paper, the evolution of reconnected tubes across a uniform CS and 
skewed magnetic fields is analyzed in detail through simulations and modeling. The inner structure 
of the GDSs consists of a thermal front where the temperature increases and most of the heating 
occurs, followed by an isothermal sub-shock where plasma is compressed. The thickness of the 
shocks, comparable to the entire length of the tubes, is calculated theoretically. 

Patchy reconnecti on seems to agree with several observations of reconnection. For example, 
Longcope et al.l (|2005l ). through EUV observations and modeling, analyzed an emerging active 
region in the vicinity of an existing one. They studied non-flaring reconnection, along a separator 
overlying the volume between them. They characterized 43 loops observed by TRACE that result 
from reconnection between the active regions, and assumed a correspondence between them and 
model flux tubes. The average diameter of the loops is 3.7 Mm with very small variation, suggesting 
a sporadic or patchy reconnection. The estimated flux of each one of them is 4 x 10 18 Mx for an 
assumed magnetic field of 37 G. Figure 19 of the paper shows an illustration of the reconnection 
model applied to the active regions where the three-dimensional CS is shown, as well as a pair of 
reconnected flux tubes created by patches of reconnection. The time history of the 171 A loops 
with very similar diameters suggests a highly intermittent reconnection process. 

What triggers reconnection is still an open question, but several promising mechanisms have 



been proposed. For example, the MHD CS tear ing mode instability (|Furthlll963l ; ISturrock 



Forbes Priest 



1983 



Priest 



1985 



Biskamplll986l ) may initiate bursty reconnection. 



1966: 



Shimizu et al 



(|2009l ) showed numerically that a two-dimensional fast magnetic reconnection in a simple one- 
dimensional CS can be unstable for three-dimensional resistive perturbations. Therefore, it is 
possible to randomly eject three-dimensional magnetic loops along a CS atop a post-flare arcade, 
resembling SADs. Catastrophe models explain how magne tic configurations are stable for a long 
time and suddenly release their energy (ICassak et al.ll2005l ). The Hall reconnection model, exten- 
sively studi e d in the Geospace Environmental Modeling(GEM) magnetic reconnection challenge 
(|Birn et al.l (|200ll ). and com panion papers), y ields energy release rates consistent with the ones 
observed in the solar corona ([Drake et al.l 120061 ). 
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Even though three-dimensional outflows can be generated from purely two-dimensional con- 
figurations, the assumption of perfectly anti-parallel pre-reconnection magnetic fields seems too 
restrictive for the complex solar corona, where field lines are likely to form finite angles across CSs. 
Further more, the add i tion o f a guide field is directly related to the temporal evolution of recon- 
nection. iDrake et al.l (|2006l ) shows that ambient guide magnetic field controls whether magnetic 
reconnection, once it is established, remains steady (anti-parallel reconnection) or becomes bursty 
(reconnection with a guide field) . Through particle simulations in a system consisting on two Harris 
CSs, they show that magnetic reconnection in a large-scale system with a guide field will break 
up into many islands, as a result of tearing instability. Guide fields with a shear angle of 127° are 
sufficient to produce secondary islands that grow to finite amplitude, producing well developed flux 
tubes. 

When the standard reconnection Petschek model is extended to include skewed fields, the orig- 
inal V-shaped "switch-off" shocks are replaced by RDs (where the field direction is changed and 
plasma is acceler ated) and slow-shocks (SSs; where kinetic and magnetic energy are partially con- 



ve rted to heat.) ( Petsche 



in 



fc Thorne 



1967 



Soward fc PriestJll982l : ISkender et al.ll2003l ): see Figure 1 



Longcope fe Bradshaw ( 201ol ) for a single flux tube as it passes through these sho cks). Therefore, 



the en ergy conversion occurs in a two-st ep process. For the reconnected tubes by lLongcope et al 



(|2009i ) and iGuidoni Longcopd (|2010l ). the energy conversion also has similar behavior: first, 
magnetic energy is converted to kinetic energy at the RDs, and later this kinetic energy is partially 
converted to thermal energy at the GDSs. GDSs are parallel shocks, closely related to SSs. 

Besides the assumption of being retracting reconnected flux tubes, there are several unanswered 
questions about SADs. To the best of our knowledge, there h as not been a satisfactory ex planation 
of why they seem to be hot and relatively devoid of plasma. iMcKenzie &: Hudson! (|1999l ) reported 
temperature ratios between the voids and the surrounding rays of the fan of 7.9 to 9.1. The wakes 
seem to have a slightly lower ratio: 8.6 to 9.1. The density ratios between the voids and the rays 
ranged between 1.3/2.3 to 1.0/4.0, depending on the assumptions about the line of sight depth. 
The darkness of the voids in X-ray and extre me-ultraviole t (EUV ) images and spectra, together 
with the lack of absorption signatures in EUV llnnes et al.l (l2003bl ) are best explained as pockets 
of very low plasma density, llnnes et al.1 (l2003bl ) ruled out the possibility that the dark tracks are 
caused by coronal rain (cold plasma falling gravitationally). 

In the present paper, through a simple but realistic fluid model of patchy reconnection, we show 
that plasma depletion naturally occurs in flux tubes that are reconnected across Syrovatskii-type 
CSs (| S vrovat skij 1 1 9 7 ll ) . In this kind of CSs, the background magnetic pressure has its maximum at 
the center of the CS plane, and decreases toward the edges of the plane. While the tube descends 
toward the top of the arcade, plasma density is decreased by tube's lateral expansions in response 
to the external background pressure. After the tubes lay on top of the arcade, rarefaction waves 
continue decreasing plasma density. 



Syrovatskii-type CSs may be ubiquitous in the solar corona. lEdmondson et al.l (|2010l ) show 
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that even a simple multipole topology would form very stable CSs bounded by two Y-type nulls, 
where an X-type null deforms under stresses due to motions in the photosphere. They did not 
include a guide field, but suggested that if added, the expected three-dimensional islands formed 
after reconnection may resemble longer flux-tube-like magnetic islands. 

There is another piece of the puzzle that has not been completely understo od: SADs 1 reported 
speeds (30 - 500 km s" 1 ) are smaller than the assumed Alfven speed in the corona ( McKenzie &: Hudson 
19991 : lMcKenziell200cL ksai et alJ[20ol : iMcKenzie & Savage! I2OO9I : ISavaee et alJboid ). and this con- 
tradicts any standa rd reconnection the o ry. A possible explanation for reconnected tube's decel- 
eration is given by iLinton &: Longcopd (j2006l ) . They showed that a drag effect from the plasma 
around the tube can reduce its speed considerably. The external fluid deforms and accelerates as 
the tube moves through it; this requires energy and therefore slows down the tube. The ratio be- 
tween the speed observed in their simulations and the perpendicular Alfven speed (expected speed 
from the TFT model) increases with the reconnection angle and ranges approximately between 0.1 
and 0.5 for the low plasma-/3 (ratio between thermal pressure and magnetic pressure) case. This 
"added mass" effect accounts for most of the flux tube drag at high plasma-/?, but is not enough, 
however, to explain the speed reduction at low plasma-/3 (the expected speeds taking into account 
added mass effects are double the one observed in the simulations). They conclude that there is an 
additional source of drag in this case, but they do not speculate about its origin. 

We do not know the correct answer to the disparity between observations and models regarding 
the speed of localized reconnection outflows, but most likely the slowing down is due to the inter- 
action between tubes and their surroundings. In our present work, for simplicity, this interaction is 
not included; the reconnected tubes are assumed to be completely isolated from their background. 
Therefore, we restrict ourselves to standard reconnection scenarios where the outflows are Alfvenic. 

The layout of this paper is as following. The first section describes the background configu- 
ration where reconnection is assumed to happen. Three different locations of reconnection across 
the CS are described. The next section presents the TFT equations that govern the evolution of 
the reconnected flux tubes, as well as analytical equilibrium solutions to the perpendicular part of 
the equations. Section H] describes simulations of the evolution of reconnected tubes for each recon- 
nection location. The parallel dynamics of the tube is described in detail before the tube arrives 
at the top of the underlying arcade. In Section 14.21 the- arrival at the arcade apex is simulated 
by an overdamped spring force exerted by the compressed arcade that slows down and stops the 
tube. After this arrival, the parallel dynamics changes and strong rarefactions develop. In Section 
14.31 the temporal evolution of the tube's total emission measure and mean temperature are shown, 
and a method is described to determine where in the CS a tube was reconnected. The last section 
presents a discussion of our results and their possible observational consequences. 
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2. Flare Current Sheet in the Solar Corona 

For the location of the reconnection episode in the solar corona, we assume a realistic magnetic 
field configuration like the one shown in Figure [2j This initial equilibrium background presents a 
CS in the x — y plane, located high in the corona; and skewed magnetic fields on eac h of its sides . 



The mathematic al description of this configuration is usually called Green-Syrovatskii (|Greenlll965l : 



Syrovatskiilll97ll ) CS, or "double y-type" due to the shape of its separatrices when projecting field 
lines in the y — z plane (see panel (b) of Figure[2]in the present paper, and Figure 3.b in Syrovatskii's 
paper). We included a constant guide field in the ignorable x direction to provide an arbitrary angle 
between the pre-reconnection field lines. We also assume uniform initial values of density p e , and 
pressure P e , on each side of the CS. The CS extends infinitely in the rc-direction, but is finite in the 
y-direction. The bottom edge of the CS marks the top of the flare arcade (shown in Figure [2] as 
arc-shaped loops) with footpoints in the solar surface (bottom panel of the box). We will denote 
the half length of the CS in the y-direction as L e . 

The entire region is assumed to be free of electrical resistivity, with exception of a short-lived 
small patch, somewhere in the CS where reconnection occurs (small sphere, in Figure [2j). The 
horizontal and vertical position of the reconnection site will be denoted by xr and yn, respectively. 
Any horizontal position of the reconnection site is equivalent due to the initial 2 1/2 dimensional 
symmetry of the magnetic configuration described above. For simplicity, we will assume it to be at 
the origin of the horizontal axis, then xr = 0. 

We will distinguish three reconnection cases, based on their location on the CS. If reconnection 
happens in the top half of the CS (0 < yR < L e ), it will be called TOP reconnection, if it happens in 
the middle of the CS (yn = 0), it will be denoted CENTER reconnection, and finally if reconnection 
occurs somewhere in the bottom half of the CS, it will be referred as BOTTOM reconnection. Figure 
[3] represents the region of the CS enclosed by a dashed rectangle in Figure [2] and shows examples 
of the three different reconnection positions. For the sake of concreteness, we have assumed that 
the half length of the CS is L e = 100 Mm , a typical length in the solar corona. 

The magnetic field near the CS can be described as the limit of z — > of the skewed Green- 
Syrovatskii magnetic field, as follows 



B e = xB ex ± yB, 



ty 1 



= SiB ex ± yB eyy /l - y>\ for z \ 0. (2-1) 

Here, B ex and B ey are constants. B ex is the guide field, and B ey is the y-component of the usual 
Green-Syrovatskii field. The positive (negative) sign corresponds to the front (back) side of the CS. 
For simplicity, we defined unitless prime- variables, y' = y/L e and x' = x/L e . The lines y' = ±1 
represent the edges of the CS. The magnitude of the external field is maximum at the center of the 
CS and decreases toward the edges, as so does the corresponding local Alfven speed. 



Only field lines in the vicinity of the plane of the CS that intersect the resistive patch will 
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be reconnected, forming two opposite thin flux tubes. In Figure O spherical resistive patches are 
represented by their circular projections in the plane of the CS. Reconnected flux tubes have a finite 
thickness, as also does the CS. The tube's cross sectional area at the reconnection region, A, is 
determined by the radius of the reconnection patch. This area and the magnitude of the magnetic 
field at the reconnection site determine the flux of the tube, 4> = BA. 

The actual size of the reconnection patches in the solar corona has not been determined, but if 
SADs represent the cross sectional area of reconnected fl ux tubes, this value can be estimated to be 



only a few megameters in radius, as shown in Figure [TJ iMcKenzie Savagd (|2009l ) reported dark 
voids sizes of approximately 10 7 km 2 , with significant range for variation. The areas of the dark 
pockets can only give an estimation because they may change as tubes move, due to pressure balance 
with the background. Typical flare coronal loops extend for a hundred megameters, therefore it is 
reasonable to assume that the reconnected tubes are thin compared with their length. 

Field lines in front of the CS that intersect a given reconnection region will connect at this 
site with field lines on the other side of the CS. Figure [3] shows a sample of three pre-reconnection 
field lines on each side of the CS that intersect a TOP reconnection site located at y' R = 3/4. We 
will assume that the pre-reconnection flux tubes are parameterized by one central field line (their 
axes, thicker lines in the graph), that crosses the center of the reconnection region. In this figure, 
only pre-reconnection representative field lines are drawn for the CENTER case. 

The half angle Cr between the pre-reconnection representative field lines at a given reconnection 
location is a parameter in our model, and can be calculated as 

tan (Ch) = ^Jl-V'i = tan (ftWl - V% (2-2) 

where Co represents the half angle that field lines make at the center of the CS. Both angles are 
shown in Figure 

Any pre-reconnection field line near the CS can be described as 

y'(x') = ±sin [tan (CoK + c] , z^O, (2-3) 

with c being a constant that can be determined by any point that intersects the curve. For field 
lines that cross the reconnection point (0,y' R ), c = ±arcsin(y^), where the top sign corresponds to 
field lines in front of the CS [z > 0, thick solid lines in Figure [3]), and the bottom sign corresponds 
to field lines at the back of the CS (z < 0, thick dotted lines in Figure [3]). 

The x-position at which a tube (its representative field line) is the closest to an edge of the CS 
can be found replacing y' = ±1 in Equation ()2-3p (the positive sign corresponds to the top edge of 
the CS, and the negative sign corresponds to bottom edge of the CS), to obtain 

\ — arcsin (y' R ) 



x' 



top 



x bottom ~F 



tan (Co) 
I + arcsin (y' R ) 
tan (Co) 



z^O, 

z^O. (2-4) 
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At these locations, the tubes are tangent to the edge of the CS since the only remaining component 
of the magnetic field is the guide field (Equation (|2-ip ) . The tubes have legs that continue downward 
to the solar surface and these sections of the tubes do not generally lay in the same plane as the 
CS (see Figure E]). 

After reconnection, the reconnected flux tubes are symmetric in the x-direction, but the initial 
2 1/2-d symmetry of the system is lost. The localization of the reconnection episode makes the 
problem three-dimensional. Fortunately, the reconnected tubes are mostly two-dimensional as they 
lay in the plane of the CS, and their retraction is along this plane (magnetic tension force is parallel 
to the plane of the CS). In Figure O representative lines for reconnected flux tubes in a BOTTOM 
reconnection site are shown. The dashed line represents the initial configuration for a reconnected 
flux tube whose left side was part of a tube on the back side of the CS and whose right side was 
part of a tube in the front side of the CS. This new tube is sharply bent near the reconnection site 
and will move upward as a transient feature that slides along the CS, between the flux layers. We 
call these tubes UPWARD moving tubes. The mixed dotted and dashed line in the figure is the 
initial configuration of the DOWNWARD moving reconnected tube that will retract in the negative 
y-direction. The directions of motion for each tube are indicated in the figure by arrows pointing 
outward from the BOTTOM reconnection site. The UPWARD and DOWNWARD moving tube 
definitions are valid for all reconnection cases. Any reconnected UPWARD moving tube is a mirror 
image of a DOWNWARD moving tube if both tubes have the same absolute value y' R reconnection 
location, but with opposite sign. They will move in opposite directions, but their dynamics and 
shapes are equal. Therefore, in this paper we will only focus on DOWNWARD moving tubes. 

The shape of a reconnected tube changes with time mainly as result of magnetic tension. This 
motion can be characterized by its representative field line arc-length parametrization, R(Z,t). At 
each point of the tube, the local tangent unit vector and curvature vector can be calculated from 
the parametrization as b = dTl/dl and k = dh/dl, respectively. Figure [3] shows an example of the 
parametrization vector for an initial UPWARD moving tube, at a generic arc-length I. Also, at this 
point of the curve, b and k are shown. Pieces of the tube like this one, that are far enough from 
the reconnection site, are initially in equilibrium since they were part of the original background 
configuration. There, the force due to the curvature of the tube is balanced by the force due to the 
external gradient of the magnetic field. 



3. Low Plasma-/3 Thin Flux Tube Equations 

Reconnection is assumed to happen in a small and short-lived region, therefore the recon- 
nected tubes are thin and their evolution can be described by the TFT equations presented by 



Guidoni k, Longcopd (|2010l ). These equations assume total pressure balance between the tubes 



and their surrounding plasma, and the tubes are assumed to be isolated from the background (no 
viscous momentum and heat exchange). The plasma in these tubes satisfies the frozen-in field 
condition since the resistivity that may have originated the reconnection episode was only non-zero 
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at the reconnection region. 



In the solar corona, the plasma is strongly magnetized, ergo its plasma-/?, is small. For 
this reason, thermal pressure can be neglected with respect to magnetic pressure. This approx- 
imation and the assumption of total pressure balance constrain the magnitude of the magnetic 
field inside the tube; at all times, it is equal to the magnitude of the external magnetic field B e 
(jGuidoni fc Longcopdl2010l ). and can be calculated from Equation 12- 1 1 



The cross sectional area of each element of the tube changes as the element moves between 
the flux layers of the CS, satisfying A = The density of each tube element of mass 5m varies 
according to p = 5mB e /(p5l, where 51 is the length of the element. 

The magnetic field direction is completely described by the positions of the tube elements 
since the frozen-in condition guarantees that fluid particles and field lines move together. These 
positions can be obtained integratin g the velocity satisfying the low plasma- f3 TFT momentum 
equation (jGuidoni Longcopdl2010l ) 

8vr 



Dv 



P 



Dt 



b 



-OP 



dl 



V 
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Here, p, v, P, and r\ are the mass density, velocity, plasma pressure, and viscosity of the fluid inside 
the tube, respectively. D/Dt is the advective derivative in the tube element's reference frame and 
Vj_ represents the gradient in the direction perpendicular to the flux tube. 

The inviscid part of the above equation has a parallel-to-the-magnetic-field term with a deriva- 
tive of the thermal pressure that cannot be neglected with respect to any magnetic pressure term 
since magnetic forces act only in the perpendicular direction (as shown in the second and third 
term of the right-hand side). The viscosity term has components in both directions but is usually 
small since viscous Reynolds numbers in the corona are very large. The exception to this rule are 
shocks, where velocity gradients become significant. For strongly magnetized plasmas, the viscosity 

5 

depends on temperature as rj = rj c T2 with r/ c being a constant. 

The contribution of thermal conductivity to the change in entropy of the plasma is approx- 
imately two orders of magnitude larger than the viscosity. The change in entropy of the plasma 
(assumed to be an ideal gas) of any given tube element can be expressed as As = In P c /P e , where 
P c is related to the tube element's pressure and density as 



P(t) = p c (t) 



Pc 



(3-2) 



Here, 7 is the adiabatic gas constant. 



The TFT entropy equation can be transformed in an equation for P c , as follows (jGuidoni &: Longcope 



2010T ) 
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The temperature dependence of the thermal conductivity is the same as for the viscosity, k = n c T^ , 
with k c being a constant. The ideal gas constitutive relation is P = pksT '/m = 2nkBT, with ks 
being the Boltzmann constant, m the average particle mass, and n the electron number density. 



The TFT mass conservation has the following expression (jGuidoni &: Longcopd 12010 ) 

D {B e \ B P ~ <9v 



One of the goals of this paper is to explain how reconnected flux tubes may become regions 
of density depletion, like in the dark voids mentioned in the Introduction. To that end, the above 
equation can be re-arranged to show that the density of each tube element may change due to three 
different effects, 

D D dvn 

-ln(p) = --ta^-gl+vx.k, (3-5) 

where v\\ is the parallel component of the velocity and vj_ is the perpendicular one. 

The first term in the right-hand side corresponds to a change in area of the tube as the tube 
element moves to regions of different confining magnetic field magnitude. In addition, if neighboring 
elements are slowing down ahead of the direction of motion, a pile-up occurs and the density of 
the element increases (its length gets reduced, and its area is fixed by the external field). The 
opposite effect occurs if the neighboring elements are slowing down upstream. The second term 
in the right-hand side reflects this effect, which is similar to what happens in traffic jams. The 
third term in the right-hand side is related to the expansion or contraction of a curved tube in its 
perpendicular direction. One example of this effect is a circular tube with fixed area expanding 
or contracting radially. In this case, the curvature vector always points toward the center of the 
tube, and the perpendicular velocity is parallel (contraction) or antiparallel (expansion) to this 
vector. Each tube element expands or contracts resulting in a change in density. We will show the 
relevance of some of these terms in section 14. 1[ 



3.1. Joined Equilibrium Solution 

A simpler magnetic background confi guration than the one assumed in the present paper 



was studied by iGuidoni Longcopd (12010 ). There, the background field is skewed and uniform, 
and the shape of the reconnected thin flux tubes, as they evolve, can be described by a joined 
equilibrium solution (JES). We will call this uniform background model GuidonilO. In this case, 
the reconnected tubes, at any instant, are composed by three straight segments joined at two 
corners (the bends) that move at the Alfven speed along the initial tube. There, the magnetic field 
is rotated without changing its magnitude (rotational discontinuity). Each straight segment of the 
tube is an equilibrium solution of the inviscid perpendicular part of the momentum TFT equation, 
but the joined solution is not in equilibrium due to the sharp angle at the bends. The two lateral 
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segments are at rest and the central segment moves at the projection of the Alfven speed in the 
y-direction. Since the magnetic field magnitude and the initial density are uniform, the Alfven 
speed is the same everywhere. 

The actual parallel dynamics in GuidonilO does not satisfy an equilibrium solution (constant 
pressure and temperature) since strong GDSs develop inside the tubes. Nevertheless, outside 
the shocks, the density and temperature are approximately constant. To completely satisfy an 
equilibrium solution of the general TFT equations, the change in velocity should satisfy b • ^ = 0; 
a condition that is not satisfied at the shocks. 

For our current non-uniform magnetic configuration (Equation (|2-ip ). it is also possible to find 
an analytical solution for the inviscid perpendicular part of the momentum TFT Equation (13-11) . 
where each mass element is in equilibrium (j£ = 0). In this case, the analytical solution differs 
quantitatively with the actual tube dynamics, as we will see in the next section, but the general 
aspects of the solution are maintained. 

It is straightforward to see that a tube satisfying the following equation constitutes such an 
equilibrium solution 



A possible general equilibrium solution of the TFT equations is given by a tube having uniform 
density and temperature that satisfies the above equation, as well as b • & = 0. In this section, we 
will only focus on the shape of the tube given by (|3-6p , as the parallel dynamics is not expected to 
be in equilibrium. 

The reconnected tubes lay in the plane of the CS, then the derivative with respect of the 
arc-length of the tube can be expressed as d/dl = (l/vl + x' 2 )d/dy' , where x' = dx' /dy'. With 
this expression, after some algebra, Equation (|3-6|) becomes 



The initial background configuration satisfies the general MHD equation's equilibrium. Its field 
lines (12-3D . also satisfy the above equation. Therefore, pre-reconnection tubes are in equilibrium 
with respect to the general TFT equations. Immediately after reconnection, every tube element is 
also in equilibrium, except the central element that is sharply bent at the reconnection site. 

The reconnected flux tubes are symmetric in the x direction, and smooth at their center. A 
solution to Equation 13-71 that satisfies these conditions can be expressed as 



k 




(3-6) 




(3-7) 





This solution is a one-parameter family of curves whose parameter is the y-position of the center 
of the tube, y' c . 
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It is possible to construct a JES where the shape of the reconnected tubes is composed of 
three different equilibrium regions: two unperturbed equilibrium side sections given by Equation 
(|2-3p . and a central curve described by the above equation. The left panel of Figure 2] shows 
such solutions for a CENTER reconnection case with half reconnection angle Cr = Co = 45°, and 
DOWNWARD moving tube. Each dotted line corresponds to a different parameter (position of 
the center of the tube, represented in the figure by triangles). The unperturbed side sections are 
part of the initial reconnected tube (solid lines in the figure). For example, a complete JES for 
the parameter y' c labeled "C" in the figure, corresponds to two side sections of the initial tube 
(sections "AB" and "DE"), plus a central symmetric curve ("BCD"). For a given center position, 
the intersection between the dotted curve and the solid line (the bends, shown as circles in the 
figure) can be obtained setting Equations (|3-8p and f|2-3[) equal to each other. 

The whole set of these JESs describes a reconnected tube that descends from the reconnection 
region. At each time, the tube is shorter than the original tube. Even though each section of the 
joined solution is an equilibrium solution, the joined solution is not. The bends are clearly out of 
equilibrium, therefore a complete curve like "ABCDE" does not satisfy equilibrium Equation (!3-7p . 
Left panels of Figures [5] and [6] show similar joined equilibrium solutions for BOTTOM and TOP 
reconnection cases, respectively. 

For comparison purposes, all cases analyzed in the present section have a half reconnection 
angle Cr = 45° (note that this is different for the CENTER reconnection case shown in Figure [3] that 
has Cr ~ 57°; this choice visually simplified the graph making pre-reconnection field lines parallel 
to each other on each side of the CS, for all the reconnection cases). From now on, CENTER 
reconnection case implies Cr = 45°. 

In the CENTER and BOTTOM reconnection cases, the concavity for the JESs is up for all 
the central regions. On the other hand, for the TOP reconnection case, the concavity of the central 
regions changes from down to up when the center position changes from positive to negative (at 
y' c = 0, the equilibrium central region corresponds to a straight horizontal line). 

Equilibrium solutions ()3-8[) are valid as long as the position of the center of the tube is above 
the edge of the CS (y' c > —1, or y c > 100 Mm). This corresponds to the arrival of the center of the 
tube to the top of the arcade. Since our analysis was restricted to field lines near the plane of the 
CS, equilibrium solutions beyond this point have no real significance. 

4. Simulations 

The low plasma-/? TFT Equations (|3-1|> . (|3-3p . and (|3-4p are non-linear, hence a time-dependent 
analytical solution of these equations is hardly a possibility. The JES, although far from being a 
general solution, presents the correct general aspects of the shape of the tube, as we will see in this 
section. The parallel dynamics are a different matter. The strong compressions of the plasma due 
to the shortening of the tube generate GDSs along the tubes that present discontinuities in the 



state variables (jGuidoni &: Longcopdl2010l ). 



To study the full dynamical solution of the TFT equations, we developed a computer program 
called Dynamical Evolution of Flux Tubes (DEFT) which solves the TFT Equations (|3-ip and (|3-3|> 
in dimensionless form for the two reconnected tubes (DOWNWARD and UPWARD moving). Mass 
conservation is automatically satisfied since the DEFT program implements a Lagrangian approach 
where each mass element of the tube is followed. The program uses a one-dimensional staggered 
mesh where each tube piece is represented by grid points at its ends. A rather detailed description 
of this computer program was presented in GuidonilO. 

In this paper, we will present three different simulations that correspond to the three recon- 
nection positions and angles described in section 13.11 The initial tube shapes correspond to a half 
reconnection angle C,r = 45°. Results depend quantitatively on the reconnection angle, but we 
chose this intermediate value to show the general aspects of the solutions that are common for a 
wide range of angles. The half angle that field lines make at the center of the CS is Co — 57° for 
the BOTTOM and TOP cases, and Co - 45° for the CENTER case. 

The central initial part of the tubes are clearly out of equilibrium due to the sharp angle at 
the reconnection site. To avoid introducing length scales at the limit of resolution, the central part 
of the initial tube is smoothed. The end points of the tube are assumed to be fixed, and to have 
no temperature gradient (this last condition ensures no heat transfer from the end points). 

For concreteness, we will present the results from all simulations in units that assume a magni- 
tude of the background magnetic field at the reconnection site of 10 G, an initial uniform electron 
number density of n e = 10 8 cm" 3 , and initial uniform temperature T e = 1 MK. The Alfven speed 
VaeR at the reconnection site is then approximately 2200 km/s. As in the above section, the half 
length of the CS is L e = 100 Mm. We assume fn equal to half the mass of the proton. Each tube 
has 1100 points. 

The unitless numbers for this simulation correspond to plasma-/? ~ 0.007, Prandtl number 
(ratio between the viscosity and thermal conductivity) P r = 0.01, and viscous Reynolds number 
R v = L e v ae Rp e /r] ~ 2982. These values are typical values for the high solar corona. 

The simulated evolution of a DOWNWARD moving tube reconnected at the CENTER site 
is shown in the right panel of Figure HI This simulation has similarities with the corresponding 
JES (left panel). For each given time, the tube consists of three sections: two unperturbed ones on 
its sides and a central one that moves downward. Times shown were chosen to coincide with the 
same bend positions as in the corresponding JES (circles are located in the same positions). The 
concavity of the central portion has the same sign as the JES, although the JESs are more curved. 

The flattening of the center part of the tube is due to the perpendicular component of the 
viscous force density in Equation (|3-ip . In the figure, the small crosses near the center of the tube 
correspond to the location of sub-shocks, where parallel velocities have large gradients and viscosity 
becomes relevant. At those locations, the perpendicular viscous force acts in the direction opposite 
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to the local curvature vector; therefore, the total force in the curvature direction is reduced. This 
increases the curvature of the tube, resulting in a change of curvature sign at the center of the tube. 
A thorough analysis of the shocks is presented in the next section. 

In the simulation, the center of the tube moves with approximately constant speed equal to 
the y-projection of the reconnection Alfven speed, which is of the order of 1500 km/s. This can 
be understood, at least for earlier times, by comparison with the GuidonilO model where this is 
also the case. The initial tube can be approximated as the initial conditions for skewed uniform 
magnetic fields (same as in GuidonilO). The magnitude of the local Alfven speed decreases toward 
the bottom edge of the CS, therefore there are sections of the tube, like the center one, whose 
perpendicular speed is super- Alfvenic. The ratio between some tube elements' perpendicular speed 
and the magnitude of the local Alfven speed can reach values above 1.1, although this value is 
relative since the direction of the background magnetic field is different than the one of the tube. 

In over a minute, the central part of the tube arrives at the lower edge of the CS which marks 
the beginning of the top of the arcade, located at y = —100 Mm. The evolution of the tube after 
its arrival at the top of the arcade is discussed in detail in Section 14.21 

In the GuidonilO model, the bends are RDs that move at the Alfven speed. We expect similar 
behavior in our present case. The positions of the circles in right panel of FigureEJwere theoretically 
calculated as moving at the local Alfven speed along the tube. They coincide with the actual bend 
positions in the simulation. They continuously slow down in the y-direction as they approach the 
edge of the CS. Their speed in the x-direction, determined by the guide field, is constant. Being 
faster than the bends in the y-direction, the tube's center is the first to arrive at the top of the 
arcade. This difference in speed is what determines the general upward concavity of the moving 
part the tube. 

At the bends, the initial stationary plasma is deflected along the bisector of the angle between 
the two adjacent tube portions (average direction of the curvature force). The resulting speed is 
Alfvenic and has a parallel component toward the center of the tube. These bends, as in GuidonilO, 
are RDs where only the magnetic field direction and the plasma velocity are changed. Notably, the 
angle at the bends remains approximately constant during the entire evolution of the tube, and 
equal to 180 — (r degrees (same as in the uniform case). 

The BOTTOM reconnection case simulation (shown in left panel of Figure [5]) presents the 
same characteristics as the CENTER one, but the time it takes for the tube to arrive at the top of 
the arcade is much shorter (~ 16 s), as the reconnection site is closer to the apex of the arcade. 

On the other hand, the TOP reconnection case (right panel, FigureEJ presents some differences. 
The general concavity sign of the tube changes with time, and this change occurs much later in the 
evolution of the tube compared to the corresponding JES (left panel). Initially, the bends accelerate 
in the vertical direction and move faster than the center of the tube, leading to a downward concavity 
(same as in JES). After the bends reach the center of the CS, they start decelerating as they move 
toward the top of the arcade. The tube's center accelerates continuously, catching up with the bends 
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after approximately 70 s, where the tube is almost flat. Meanwhile, the perpendicular viscosity 
force density increases the curvature near the locations of the shocks, curving the center part in the 
opposite direction with respect to the general curvature of the tube. For this TOP case, the center 
of the tube is not the first part to arrive at the top of the arcade. The ratio between some tube 
elements' perpendicular speed and the magnitude of the local Alfven speed reaches values above 
1.7. 



4.1. Parallel dynamics 

The parallel dynamics along the reconnected tubes presents complex features like shocks and 
coupled non-linear waves. There is no simple equilibrium counter-part to compare it with, as there 
was for the perpendicular motion. Nevertheless, its complexity is what makes the parallel dynamics 
very rich in behavior, and this section is devoted to its description. 

The bends accelerate plasma supersonically toward the center of the tu be where they collide 



and ar e stopped abruptly. As a result, two strong GDSs are launched outward. iGuidoni fc Longcope 



(|2010i ) describe in detail their inner structure for the uniform and skewed magnetic field config- 
uration. It consists of one thermal front, where most of the heating occurs, and an isothermal 
sub-shock, where most of the compression occurs. The thickness of the GDSs increases with time 
and can be of order of tens of megameters, which is comparable to the size of the entire reconnected 
tube and larger than the particle's mean free path. The large and small plus signs in the right 
panels of Figures 01 [5j and [6] represent the beginning (heat front) and the end (sub-shock) of the 
GDSs, respectively, for each simulation time. 

For the three reconnection cases presented in this paper, the GDSs are ahead of the bends at 
early times after reconnection. Initially, the shocks develop their inner structure quite rapidly until 
they achieve a near-equilibrium thickness. Meanwhile, the central part of the tube is increasing 
in length due to the difference in vertical speed between the center and the bends. Therefore, the 
central portion eventually becomes larger than the distance covered by the shocks, and the bends 
move ahead. For the TOP reconnection case, this situation is reversed after the bends move to the 
lower half of the CS, where they slow down and the shocks overrun the bends again. 

The electron density has a sharp increase at the sub-shocks where plasma is extremely com- 
pressed. Here, the second term of the right hand side of Equation (|3-5p is the most relevant to 
determine the electron density. The top panel of Figure [7] shows density profiles for the CENTER 
reconnection case where each solid line corresponds to a different time in the simulation. At early 
times, the density at the center of the tube is over six times the background density. Even though 
the theoretical density jump limit for steady state GDSs is 4, we achieve higher values in our sim- 
ulations because there is an initial transient over-shoot of density at the center of the tube. The 
plasma requires a finite time to adjust to the correct steady state shock strength. Density decreases 
as times goes by, although never reaching the equilibrium value. In the figure, we include an inset 
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to the left that reproduces the t = 60 s electron density profile, where each point of the graph 
corresponds to a grid point of the tube. This shows that shocks are well-resolved by the DEFT 
program. 

Density enhancements by shocks are a common feature in any standard reconnection model; 
this compression occurs at slow-mo de shocks or switch-off shock s, which are close relatives of the 



gas-dynamic shocks of our paper (jLongcope fc Bradshawl |2010| ). Therefore, dense outflows are 



usually expected from reconnection. In our patchy reconnection case, central long-lasting, high- 
density post-shock regions, or plugs, expand and move outwardly from the reconnection site. Their 
increase in size can be seen in the right panels of Figures [U [5l and [6J as tubes descend. The 
section of the tubes limited by the small crosses corresponds to the plug. This feature should be 
observable by an X-ray instrument of sensibility suitab le to the emission measure. In fact, bright X - 



ray features apparently shrinking have been observed ( McKenzid 12000 : iMcKenzie &: Hudso 



origntA- 

Aooih . 



Sheelev fc Wand (|2002j ) also describe some coronal inflows as barely visible density enhancements 



that move sunward, with dark tails forming behind them. 

Between the bends and the sub-shock, the electron density decreases before increasing again 
at the sub-shocks. This decrease in density is caused by the lateral expansion of the tube as it 
moves into regions of lower magnetic field (first term in the right-hand side of Equation (|3-5j) ). as 
well as by the expansion of the curved tube in its perpendicular direction, although the latter effect 
is smaller than the former. In the figure, the positive ^-positions of the bends are shown as circles 
for each chosen time. After 60 s, the plasma density in these regions is less than 70 percent that of 
the background density. If observed, this tube would present dark and elongated thin regions (low 
emission measure) that get darker, descend toward the arcade, and grow in length as the bends 
move along the legs. This resembles the SADs or sinking columns mentioned in the introduction. 
By the time the tube arrives at the top of the arcade, it would be mostly dark, and approximately 
half of the dark region would be cold, and the other half would have a range of higher temperatures 
as the heat fronts move outwardly. 

The GDSs are strong and the temperature jump across them may exceed an order of magnitude. 
The CENTER reconnection case's temperature profiles are shown in the middle panel of Figure [7| 
for the same times as in the top panel. The temperature across the shocks initially increases with 
time, achieving values higher than 7 MK, and later decreases. The large and small cross signs are 
the positions of the beginning and end of the thickness of the shocks, respectively. The length of the 
thermal fronts increases with time as they move outwardly from the center of the tube. The inset 
to the left reproduces the t = 60 s temperature profile, where each point of the graph corresponds 
to a simulation point in the tube. 

The thermal pressure gradient in Equation (|3- 1 1) is crucial for the generation of the shocks. If 
it were erroneously neglected in this equation, the plasma would continuously pile up at the center 
of the tube without limit. However, this term is not the only factor in the determination of parallel 
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dynamics. The parallel, inviscid part of Equation ()3-l j) can be re-written in the following way 

Dv\\ Db ldP / . \ 

~Dt = v ^^-p5T' (4 " 1} 

showing two possible sources for the change in tube element's parallel velocity (except perhaps 
inside the sub-shocks where the viscosity contribution may also be large). The first term on the 
right-hand side is a fictitious force caused by the change in direction of the unit parallel vector as 
tube elements move. The solid dark squares in Figures [U and [6] show the position of a generic 
tube element as it moves along the tube. 

The parallel velocity profile for the CENTER reconnection case is shown in the bottom panel 
of Figure El for the same times as in the other two panels. The positive x-positions of the bends 
are also shown as circles for each chosen time. There, the plasma at rest is abruptly accelerated 
to super-sonic speeds, therefore the slope of the parallel velocity curve is almost vertical at these 
points. The same occurs at the bends on the other side of the tube where the inflows are directed 
to the right (positive speed). Since the angle at the bends remains mostly the same as time goes by, 
the value of the speed achieved at the bends is also approximately constant. The last time shown 
in this figure is 60 s, before any part of the tube has arrived to the top of the arcade. 

The parallel velocity after the bends slightly increases due to a small predominance of the 
first term in Equation (|4-ip . This term also non-linear ly couples the parallel dynamics (sound 
waves) to the perpendicular one (Alfven waves), which introduces dispersion. The oscillations in 
parallel velocity that follow the bends are not an artifact of the simulation; they are real dispersive 
waves due to this coupling. This effect was not present in GuidonilO since the tube segments 
were straight at all times. These waves are smoothed out by diffusive processes, mostly related to 
thermal conduction. 

The following slow parallel deceleration coincides with the thermal front positions, where 
the second term of Equation (|4-ip becomes relevant. There, the inflows are decelerated when 
tube elements encounter the heat front, until they are suddenly stopped at the sub-shocks (sharp 
decrease in parallel velocity near the center of the tube). In the post-shock region, the flows are 
reversed in direction with respect to the post-bend velocity. 

The BOTTOM reconnection case presents the same general features as the CENTER recon- 
nection case for the three quantities plotted in Figure [71 These two cases are similar because their 
dynamics corresponds purely to the lower part of the CS where the background magnetic field 
decreases monotonously toward the edge of the CS. 

The TOP reconnection case presents some differences with respect to the other cases. The 
three panels of Figure [8] show the same type of graphs as in Figure [TJ for the TOP reconnection 
case. The post-shock electron density (almost one order of magnitude higher than the ambient 
density) increases while the plug is in the top part of the CS (the tube is being squeezed by the 
ambient pressure) and decreases after crossing the center of the CS. The squeezing increases the 
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electron density after the bends while those regions are in the top half-plane of the CS, and a 
pronounced decrease in density occurs after the bends cross the center of the CS. Those regions 
below the center of the CS expand as the ambient pressure decreases. The decrease in density after 
the bends can be half that of the background density by the time the tube arrives at the arcade. 
As the tube descends, a very dense and hot central region develops and moves downward toward 
the top of the arcade, while very dark and extended regions grow in length to its sides. 

The maximum post-shock temperature is high, over 16 MK. The thermal fronts are almost as 
long as the entire tube, as shown in the middle panel of the figure. There, the fictitious force is 
several times larger than the pressure gradient force, and has opposite sign. Plasma is continuously 
being accelerated toward the center of the tube after the bends, and stopped at the sub-shock. For 
a given tube element, it is possible to roughly estimate the value of the fictitious force, for a short 
period of time, as (|vj_|) |A0|, where (|vj_|) represents the average magnitude of the perpendicular 
velocity in the considered period of time, and A8 is the corresponding angle change in the parallel 
unit vector. We have estimated this term for this case, and it is several times larger than the 
corresponding pressure gradient term in Equation (|4-ip . Here, not only the perpendicular velocity 
is higher, but also the change in angle of the unit vector as tube elements move is larger. 

The speed achieved at the bends is almost unchanged over the entire simulation as the angle 
at the bends does not change considerably through the simulation. We have run simulations with 
different reconnection angles and this angle seems to always remain approximately equal to 180 — (j? 
degrees. The post-shock parallel velocities are reversed with respect to the direction of the post- 
bend velocities, and the resulting speeds are high, almost equal in magnitude with respect to the 
speeds achieved at the bends. 

A zoom near the bends in the right panels of Figures [U [SJ and [6] would show small oscillations 
in the shape of the tubes. There, the fictitious term is much larger than the second term in the 
right-hand side of Equation (|4-lj) , as the RDs do not change the pressure of the plasma. 



4.2. Top of the Arcade 



In the previous sections, the analysis of the evolution of the tubes was done for times when 
the tubes were above the edge of the CS. The DOWNWARD moving tubes will arrive at the top 
of the arcade after a finite time. For example, the middle point of the CENTER reconnection 
tube arrives at the top of the arcade shortly after a minute of travel. This a rrival is usually 
referred in the literature as the place of the fa st-mode shock or termination shock (jForbes fc Priest 
19831 : iForbes Malherbdll986l : lForbeslll986l ). Here, the reconnected flux tubes collide with the 
underlaying arcade. 

It is likely that the parallel dynamics will continue evolving along the legs of the tube after the 
tubes had arrived to the top of the arcade. The legs extend downward to the solar surface and do 
not lay in the same plane of the CS, therefore the simple two-dimensional analysis we have done so 
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far for the evolution of the tube cannot be applied. An exact analysis of the dynamics of the tube 
after the tube arrives at the edge of the CS escapes the scope of this paper. Nevertheless, we can 
make some approximations to study it. The decrease in length of the tubes after they arrive at the 
top of the arcade is probably much smaller than the decrease they undergo since reconnection, as 
the underlying arcade prevents them from shortening. By the time the tubes have arrived to the 
top of the arcade, their total lengths decreased by 7, 15, and 28 % for the BOTTOM, CENTER, 
and TOP reconnection cases, respectively. They may also change their length, only by a small 
fraction, by decreasing the angle between the legs of the tube and its flat central part (panel (a) of 
Figure shows a cartoon where a reconnected tube (thick line) is laying on top of the underlaying 
arcade). 

To simplify the analysis, we can assume that this length is fixed after the arrival at the arcade 
apex. The perpendicular motion of the tube is assumed to have stopped, and only the parallel 
dynamics along the legs would continue, mostly independently of how curved the tube is. To study 
this parallel dynamics, we will approximate the tube legs as straight line extensions on the sides of 
the tubes at the edge of the CS and parallel to the CS plane. Therefore, in our simulations an extra 
straight segment is added to the end of the tubes to simulate the legs; then the tubes are extended 
in the horizontal axis near the edge of the CS, as shown in the left panel of Figure H] (dashed lines). 

Most of the ob s erved dark vo i ds slo w to a stop when they arrive at the top of the arcade 



(jSheeley et al.ll2004l ). iLinton et al.1 (J2009) show through three-dimensional MHD simulations that 
reconnected flux tubes decelerate rapidly when they hit the Y-lines of a skewed Sirovatskii CS and 
the sheared arcade beyond them. Therefore, we chose in this first analysis to ignore the possibility 
of full "bounce" back of the tube toward the reconnection site. We propose a scenario where tubes 
are slowed down and stopped there by a damping force. We simulate this effect as a restoring force, 
as shown in panel (c) of Figure [H The underlying arcade gets compressed and lightly curved due 
to the arrival of the reconnected tube, and these stronger field lines exert an upward force on the 
tube. We add a perpendicular damping and spring force densities to Equation (|3-ip . active only 
below the edge of the CS 

f = pn(n 2 y-un- v) for y' < -1. (4-2) 

Here, n indicates a unit vector perpendicular to the tube, J7 2 is the spring constant, and v is the 
damping coefficient. This force only acts in the perpendicular direction, slowing down the tube and 
stopping it. 

The above force density will result in overdamped oscillations when the ratio between the 
damping coefficient and the spring force, £ = v/20,, is larger than one. We opted for this in our 
simulations in Section HI The cases £ = 1 and £ < 1 correspond to critical and under-damped 
oscillations, respectively. 



Different arcade damping coefficients lead to different damping times, but do not considerably 
affect the parallel dynamics of the tube. We ran simulations with coefficients differing by orders 
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of magnitude and all results were similar. For the results presented in this paper, the simulations 
were carried out with coefficients equal to £ = 2.0 and SI ~ 0.49 s . 

In our simulations, the reconnected tubes retract from the reconnection site to the top of the 
arcade, and slow down due to the damped spring force from the arcade, until they lay flat on top 
of it. Then, the final shape of a DOWNWARD moving tube would be a straight line on top of the 
arcade (see right panels of Figures [U and [6] for late times in the simulation) . 

Since the bends decelerate in the lower part of the CS, they are not the first sections of the 
tube to arrive at the top of the arcade. The first regions to do so become flat and grow in size 
as more tube elements descend atop the arcade (compare, for example, t = 65 and 70 s of the 
CENTER simulation in Figure H]). We will call the outer edges of these flat regions "second bend". 

Some of the parallel velocity change is due to the perpendicular dynamics captured in the first 
term in the right-hand side of Equation (|4-ip . After a section of the tube arrives at the top of the 
arcade, its perpendicular velocity is suddenly decreased by the damping spring force. Therefore, the 
first term of Equation (14- ID becomes unimportant compared to the other terms in the equation. The 
top panel of Figure [10] shows parallel velocity profiles for late times in the CENTER reconnection 
simulation (only the positive x-axis side). The dashed line corresponds to simulation time equal 
to 60 s -the earliest time shown in this graph and the last time shown in Figure [7] just before any 
part of the tube arrives at the arcade. The sharp discontinuity on the left shows the deceleration 
at the sub-shock. The next times shown in the graph present a new discontinuity at the locations 
of the second bend, where the parallel velocity is sharply decreased by the sudden reduction in 
perpendicular speed. The "SB" arrows point toward the second bend locations at t = 65 s and 
t = 70 s. 

When the second bend and the original bend become close (see arrow labeled "SB-B" in Figure 
[TO]) , the angle at the latter one decreases until it finally becomes 180° when the bend arrives at 
the top of the arcade. Then, the inflow speed is drastically reduced because it depends strongly 
on this angle. The gas becomes rarefied and the electron density decreases rapidly behind the 
bends. Here, previously accelerated plasma continues moving toward the center of the tube leaving 
a density depletion behind. By the last time of our simulation, the electron density decrease is 
more than 40 percent of the background density. The bottom panel of Figure PTOl shows the electron 
density evolution for the same times shown in the top panel. The almost vertical slopes to the 
left indicate the sub-shock positions that continue moving outwardly, even after the bends arrive 
at the top of the arcade. For the TOP reconnection case, the maximum electron density decrease 
achieved by our simulations is more than 50 percent of the background. 

The sub-shocks will continue moving outwardly and interact with the vacated regions. Al- 
though our simulations stop before this occurs, we hypothesize that the electron density will con- 
tinue decreasing for \x\ < \xb u m\ an d the central hot plug will mix with these rarefaction waves, 
generating secondary rarefaction waves that would move toward the center of the tube and ulti- 



mately will disassemble the hot plug, as described by lLongcope et al.l (j201dl ) . 
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4.3. Emission Measure and Mean Temperature 

The differences in dynamics between each reconnection case suggest an observational signature 
of where in the CS a flux tube had been reconnected. The tube's average temperature, weighted 
by the differential emission measure per flux, DEM(T), can be calculated as 

J2 DEM(T) T AT 

^ = £ DEM(T) AT ' ^ 

T>T e 

where the DEM (T) is computed for small temperature bins of size AT along the tube, as 

\T(l)-T\<^f 

DEM{T) = . (4-4) 

Only tube elements whose temperature is inside the range of the bin are included in the sum. (T) 
refers to the heated material of the tube; it excludes segments at or below ambient temperature. 
Figure [TT] shows the time evolution of the tube's average temperature (dotted line) for the CENTER 
reconnection case. 

The emission measure per flux of a given tube element, EM(1), is computed in the following 

way 

SM( . ,4-5, 

With the above definition, the tube's total emission measure per flux becomes 

EM 



T{l)>T e 

The above quantity considers only the emission measure of the heated parts of the tube (T > T e ), 
excluding cool sections like the ones following the bends and preceding the heat fronts, shown in 
the middle panel of Figure [7j Ambient plasma emission measure is not included either. Figure 
[TT1 shows the time evolution of the tube's total emission measure per flux (dashed line) for the 
CENTER reconnection case. 

In the figure, the arrow labeled "C" indicates the time at which the center of the tube arrives 
at the top of the arcade, and the one labeled "B" represents the arrival of the bends at the top of 
the arcade. Initially, the mean temperature of the tube increases at a faster rate than the emission 
measure. They both achieve a local maximum before time "C". The following negative slope for 
both curves indicates that these quantities would rapidly decrease if the tube were left to retract 
without the arcade interrupting its motion. Remarkably, they increase after this point. The length 
of the tube does not change considerably while laying on top of the arcade, but the heat fronts 
continue moving along the tube, increasing its average temperature. Our simulations end before 
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the heat fronts arrive at the footpoints, but we hypothesize that the average temperature of the 
tube will continue increasing until they do so. 

The central h ot plug lasts longer th an a conductive cooling time of a region with similar tem- 
perature and size (ILongcope et al.ll2010l ) as it is maintained by the inflows. These quantities are 
therefore long-lasting. The TOP reconnection case presents some differences. The mean tempera- 
ture and emission temperature grow together, as shown in Figure PT2l Here, the arcade also prevents 
a decrease in both quantities. 

In order to provide a general method to determine where in the CS reconnection happened for 
a given flux tube, we normalized the mean temperature and emission measure by their maximum 
achieved value. Figure [13] shows the dependence of normalized emission measure with normalized 
mean temperature for each of the reconnection cases. The BOTTOM and CENTER case (dotted 
and dashed lines, respectively) are very similar to each other. For these cases, the mean temperature 
achieves its maximum before the emission measure does. On the other hand, the TOP case presents 
a different dependence. Both quantities initially increase simultaneously and curl at the top end of 
the curve. 



Discussion 



We have presented a model of tube dynamics, accompanied by simulations, following transient 
and localized magnetic reconnection in a realistic coronal background configuration (Sirovatskii 
CS). We have shown, that retracting reconnected flux tubes may present elongated regions devoid 
of plasma, as well as long lasting, dense central hot regions. The latter are created by GDSs at the 
center of the tube, consisting in long thermal fronts followed by an isothermal sub-shock. In general, 
the jump in density across the shock exceeds the maximum value predicted by Rankine-Hugoniot 
(jRankind ll87Q ; iHugoniotl Il887l ) conditions. These jump conditions are calculated assuming the 
shock is in steady state, which is not the case here. For the TOP reconnection case, the jump in 
density can be almost an order of magnitude. This descending plasma plug, although very thin, 
would be extremely bright. It is long lived compared to a free expansion of similar temperature 
and size region since it is maintained by the inflows generated at the bends. 

GDSs are also present in reconnected tubes sliding through CS with uniform skewed fields 
(jGuidoni fc Longcopejl2010l ). However, in our present model, these hot plugs respond to the change 
in background magnetic pressure that compresses or expands them laterally, depending on which 
region of the CS the tube is sliding through. For Sirovatskii CSs, the magnetic field decreases 
toward the edges of the CS (Y-points). There, its magnitude is minimum (it is zero if there is no 
guide field). Therefore, it is possible for the tubes to move through regions of decreasing magnetic 
field as they descend toward the Sun. This seems counter-intuitive because usually magnetic fields 
in the corona are assumed to increase at lower heights. 



We have presented only three reconnection locations with a reconnection half angle of 45 °, for 
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illustrative effects. We used generic dimensional values, but results scale with B e , p e , and L e . For 
all the cases presented, the maximum plasma-/? achieved never exceeded unity. 

RDs (bends) move at the local Alfveen speed along the legs of the tubes. Elongated low density 
regions are generated by lateral expansions behind the bends, when tubes move through the lower 
half plane of the CS. There, background magnetic pressure decreases toward the edge of the CS. 
The amount of density depletion depends on how curved the tube is, which is directly related to the 
perpendicular gradient in the background magnetic field (reconnection angle). On the contrary, in 
uniform background fields, as shown in GuidonilO, tubes present only straight sections, regardless 
of the reconnection angle. No density depletion is expected in this case. When there is no significant 
density depletion, and the GDSs are allowed to extend for a detectable length, the only observed 
signature of these tubes would be a descending hot plug with high density. 

The achieved decrease in density could be as much as 30 % to 50 % of background (pre-flare) val- 



ues. This level of depletion agrees with observations (jMcKenzie fc Hudson! 1 1999 ; ISheelev k, Wang 



2002). These percentages are lower bounds since we have not considered any increase in the sur- 
rounding plasma density by, for example, chromospheric evaporation from previous reconnection 
episodes. 

The appearance of the retracting tubes changes considerably with the direction of the line 
of sight and with the location of the reconnection episode. If reconnected tubes (their observable 
parts) are seen from the same view as in panel (a) of Figure [2j they would have loop-shapes, 
although their orientation would be perpendicular to the arcade since they move in a plane that 
is parallel to the arcade axis. On the other hand, if the reconnected tubes are seen from the view 
shown in panel (b) of Figure [2j they would be seen as "hairpins" because they are tangent to the 
plane of the CS. 

For the TOP reconnection case, the concavity of descending loops when they move in the top 
part of the CS plane is the one that would be expected from a cusp-shaped loop. In this region, 
these loops would be bright at their sides (they are being squeezed by ambient plasma) and even 
brighter at their center. If the line of sight is parallel to the x-direction, this loop would appear like 
a bright small region, descending at Alfvenic speeds. After the bends cross the center of the CS, 
dark regions will develop behind them, and this view will change to vertical dark regions preceded 
by a bright region, until the center catches up with the bends and the situation is reversed. The 
dark and bright regions are hot since the heat fronts move at speeds similar to the bend speeds. 
Their sizes are comparable to the entire tube. 

For BOTTOM and CENTER reconnection cases, their concavity is always U-shaped. If seen 
from the x direction, a bright region is followed by an elongated vertical dark one that might be 
hot or not, depending on the location of the heat front. 

Ob served coronal inflows can be bright or dark, or both; they may a lso have loop- or tadpole- 



shapes (jMcKenziell2000l : lMcKenzie & Hudsonll200ll : ISheeley k, Wa ng 2002). It is possible that many 



of them are the manifestation of the same three-dimensional phenomena observed from different 
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lines of sight. The location of the tube with respect to the CS also matters. If an observation occurs, 
for example, while a tube is in the top half-plane of the CS, only bright regions would be observed. 
On the other hand, if a tube is moving through the lower portion of the CS, it would present 
dark regions, sometimes preceded by bright regions, or sometimes followed by bright regions, as 
described above. 

A dense, hot region may not always been observed due to its small size, or its high temperature. 
For reconnection angles larger than 45°, it is possible to achieve post-shock temperatures that are 
much higher than the usual instrument passbands. On the other hand, after the tube's arrival at 
the top of the arcade, this region continues growing in time, and may become visible. The pile up 
of several newly reconnected tubes may also increase the size of this region to observable size, and 
would manifest as a stationary bright region on top of the arcade. 

Remarkably, signatures of reconnection persist in a single tube, longer than the Alfvenic transit- 
time required for the tube to relax. As the tubes move downward, they encounter the top of the 
arcade that lays beneath the bottom edge of the CS. Here, the downward motion is halted by 
the arcade. We simulated this effect with a perpendicular damped spring force exerted by the 
compressed arcade. With this force, the tube comes to rest on top of the arcade as a straight line. 
The parallel dynamics continues along the legs. After they stop, the tubes cannot decrease their 
length any longer, and the RDs are shut down. Then, the gas gets rarified even more because the 
already accelerated plasma continues moving toward the center of the tube, increasing the density 
depletion behind it. 

Although our simulations stop before this, we hypothesize that the hot post-shock regions will 
continue moving along the tube and interact with the rarefaction waves, subsequently disassembling 
the hot plug. The legs of the tube would also brighten up as the thermal fronts descend toward 
the footpoints. 

The tube's parallel and perpendicular dynamics are non-linearly coupled to each other, which 
can be seen, for example, in the oscillations of the parallel velocity profiles. In addition, when the 
tube arrives at the arcade, the perpendicular dynamics is halted, and changes in the parallel speed 
are evident at the "second bend". The parallel velocity profiles are determined by pressure and 
fictitious forces due to the motion of tube elements along curved paths. This can have important 
consequences for Doppler-shifts observations. For example, the post-shock region in the TOP 
reconnection case has very strong flows in a direction almost parallel to the surface of the Sun. 
This would manifest as coincident blue- and red-shifts along the line of sight, if seen from the 
x-direction, and perpendicular to the line of sight if seen from the z-direction. Perpendicular 
velocities can be higher or lower than the local background Alfven speed, which can be related to 
the generation of secondary shocks as the tubes move (the surrounding plasma may be shocked by 
the passing tubes). 

We described the temporal behavior of the total emission measure and mean temperature of 
the heated parts of the tube and provide an observational method that may indicate where in the 
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CS the tube has been reconnected. Tubes that were reconnected in the TOP half of the CS present 
total emission measure that grows simultaneously with the mean temperature of the tube. On the 
other hand, tubes that have been reconnected in the BOTTOM part of the CS achieve a maximum 
in their mean temperature much earlier than the emission measure does. For all the reconnection 
cases, the emission measure and mean temperature duration is extended by the arrival of the tubes 
at the top of the arcade. After the arrival, the length of the tube remains the same but the heat 
fronts continue heating the plasma as they move along the legs. 

If observed downflows (SADs and SADLs) are reconnection outflows, the question of why their 
speeds are lower than the assumed Alfven speed still remains. We restricted ourselves to standard 
reconnection scenarios where the outflows are Alfvenic. Nevertheless, if dark voids are related to 
the dark regions generated behind the bends, their speeds near the edge of the current sheet are 
decreased considerably from the one they had at the center of the CS. For instance, in the CENTER 
reconnection case, the bend speeds are reduced by a fourth of what it was at the reconnection site. 
This decrease in speed of the bends is not enough to explain the observed descending voids that 
move at half of the presumed reconnection Alfven speed. However, if the half reconnection angle 
is greater than 60°, the ratio of the Alfven speed at the center of the CS to the edge of the CS is 
larger than 1/2. 

Hard X-ray sources have been ob served near top of SXT arcades, at the same time as the 



SADs descend to the top of the arcade (lAsai et al.ll2004l ). This may be the result of the interaction 



between downflows and the fast shock, or could be related to reconnection at the arcade apex 
between the retracting tubes and the arcade. Apex reconnection is a possibility in our three- 
dimensional model, but one we have not yet explored. Field lines that arrive at the top of the 
arcade have a different angle than the underlying arcade (panel (a) of Figure [9]) , and arrive there at 
Alfvenic or super-Alfvenic speeds. This is not an option in purely two-dimensional models where 
the new reconnected field lines arrive at the arcade with magnetic field direction parallel to the 
arcade. 

Maintaining a fast-shock on top of the arcade may not be feasible at sub- Alfvenic speeds. The 
standing fast s hocks exist a s long as the reconnection jets exists, therefore it is a feature of a flare 



gradual phase (|Forbeslll986l ). unless reconnec tion continues du ring the decay phase. Although some 



SADs have been reported during this phase (jAsai et al.l 12004 ). most of them were seen during the 



decay phase. It is not clear if the termination shocks would still be located at the top of the SXT 
loops if only SADs are the only reconnection outflows. 

Another possibility could be that the termination shock is located further up from the top of 
the SXT loops. If there are much hotter loops above the SXT loops, that cannot be seen in pass 
bands at lower temperatures, the termination shock could be located higher in the corona. In this 
case, the downflows could be localized downstream flows from the fast-shock, which would explain 
their low speeds. The models for fast shocks are two-dimensional and it is not clear how this would 
apply to localized and time-dependent flux tubes like the ones described here. The down side to 
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this explanation is that a fast shock certainly would compress the plasma downstream instead of 
vacating it. 

We have not consider other interactions between the reconnected tubes and their ambient 
plasma, besides pressure balance between them. The addition of a drag force to our equations 
could be an important improvement of our model that would contribute to the decrease in speed, 
as well as to make the model more realistic. We expect to include this effect in future work. The 
assumption of initial uniform density could also be improved by assuming a stratified atmosphere. 
Here, tubes that have been reconnected in the cor ona will descend in regions of higher density, and 
seen darker than the background, as suggested bv lSavage et al.l (|2010l ). although this effect has not 
been corroborated by theoretical analysis. 

We have neglected any non-fluid effects and therefore neglected non-thermal particles and their 
acceleration. 

We do not claim that the depleted regions from our model can completely explain the observed 
dark v oids, but it is suggestive tha t several observed phenomena have some similarities with our 
model. iHudson McKenzid (120011 ) stated that "Dark outflows seem inconsistent with the idea of 
heating by reconnection" . We believe our work demonstrates the contrary, and thus makes the role 
of reconnection in flares still more plausible. 
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Fig. 1. — : April 21, 2002 flare. Left panel: rotated TRACE 195 A image, taken at a time where dark voids 
can be seen descending toward the flare arcade. The dashed rectangle encloses a handful of them. Right 
panel: time sequence for the window delineated by the white dashed rectangle in the left panel. The white 
arrow in the 1:43:33 snapshot shows the position of one of the dark pockets. The other arrow points to a 
second localized void. 
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(b) 



Fig. 2. — : Flare current sheet geometry. Panel (a): field line configuration for a Green-Syrovatskii current 
sheet with a guide field in the positive x-direction. The transparent gray rectangle represents the current 
sheet in the x— y plane, located above the post-flare arcade (gray parabolic loops at the bottom). The current 
sheet is finite in the y-direction, but extends infinitely in the ^-direction (2 1/2 dimensional symmetry). The 
bottom plane of the box represents the solar surface. Black lines depict some magnetic field lines on the 
side of the current sheet that is closer to the viewer (z > 0), and the gray ones whose endpoints touch the 
top side of the box correspond to field lines on the back side of the current sheet (z < 0). The small sphere 
in the current sheet shows a generic patchy reconnection region. Field lines that intersect this region form 
a small bundle that reconnects. The dashed rectangle in the plane of the current sheet corresponds to the 
region shown in Figure EH Panel (b) : different view of the same field lines in Panel (a) , showing the double 
Y-type configuration. 
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Fig. 3. — : Excerpt of the plane of the current sheet, with three possible reconnection sites. This graph 
corresponds to the dashed box in Figure [2j Pre-reconnection field lines in front of the current sheet (z > 0) 
are represented by solid lines, and pre-reconnection field lines on the back side of the current sheet (z < 0) 
are shown as dotted lines. The edges of the current sheet are located at y = ±L e = ±100 Mm. The top circle 
shows a TOP reconnection site at yn — jL e . The set of field lines that intersect this region (pre-reconnection 
tube) is represented by a central line (thicker). At this reconnection site, representative field lines on each 
side of the current sheet are skewed by 90°. The half reconnection angle (£r = 45°) is shown to the side of 
the TOP reconnection circle. The center circle represents a CENTER reconnection site with ~ 57°. For 
this site, only the pre-reconnection representative field lines are drawn. The circle in the bottom half of the 
CS plane, located at ya = — jL e , represents a BOTTOM reconnection site with C,r = 45°. Here, the initial 
shape of the already reconnected tubes are drawn, and arrows indicate their direction of motion. The dotted 
and dashed line shows a reconnected tube that will move downward (DOWNWARD moving tube) due to 
magnetic tension at its center, and the dashed reconnected tube will move upward (UPWARD moving tube). 
For this last tube, a parametrization vector R(Z, 0) is shown for a given arc-length I. Also at this arc-length, 
the unit vector b is shown, as well as the curvature vector k. For all three reconnection cases, Co — 57°, 
therefore field lines on each side of the current sheet are parallel to each other for all reconnection positions. 
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Fig. 4. — : DOWNWARD moving tube for the CENTER reconnection case and reconnection angle (r = 
(o = 45° (note that j/-axis and cc-axis are not drawn to scale). Both panels share the same y-axis. Left panel: 
solid line represents the shape of the reconnected tube at t = 0. Each dotted line corresponds to a different 
JES with center position parameter y' c = 0.0, -0.08, -0.15, -0.23, -0.30, -0.38, -0.47, -0.55, -0.64, -0.79, 
and —0.98 (shown as triangles), respectively. Circles indicate the bends for the right side (z < 0) of the tube. 
Letters "ABCDE" join a complete JES curve for y' c = —0.79. The horizontal dashed lines show extensions 
of the tube that simulate the legs that connect it to its footpoints at the photosphere. Right panel: DEFT 
simulation of the same initial reconnected tube as in left panel. Each central curved line represents the 
shape of the tube at a different time. The time interval between each curve is approximately five seconds. 
Circles represent the theoretical positions of the bends (assumed to move at the local Alfven speed). Time 
in seconds is shown for some selected tube configurations to guide the eye. Large plus signs at the center of 
the tube indicate the beginning of the gas-dynamics shocks (heat front) and the smaller plus signs indicate 
the end of the shock (sub-shock). The solid black squares indicate the positions of a tube element as it 
moves along the tube. The arrow points to the location where the tube becomes tangent to the edge of the 
current sheet (xbottom) for the front side of the tube (z > 0). To the left of this point, the tube is extended 
as a straight line to simulate the leg that connects it to the photosphere (it corresponds to the horizontal 
dashed line to the left in the left panel). A similar straight line is added to the other end of the tube. 
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Fig. 5. — : DOWNWARD moving tube for the BOTTOM reconnection case and reconnection angle £fl = 
45° and Co — 57°, with the same format as Figure SJ Left panel: the center position parameters are 
y c l = -0.75, -0.78, -0.81, -0.84, -0.87, -0.90, -0.93, -0.97, and -1.01, respectively. Right panel: the time 
interval between each curve is approximately two seconds. 
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Fig. 6. — : DOWNWARD moving tube for the TOP reconnection case and reconnection angle Cr = 45° 
and Co — 57°, with the same format as Figure |4] Left panel: the center position parameters are y' c = 
0.75, 0.67, 0.59, 0.50, 0.40, 0.29, 0.17, 0.05, -0.08, -0.27, -0.43,, and -0.67, respectively. Right panel: the time 
interval between each curve is approximately five seconds. 
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Fig. 7. — : DOWNWARD moving tube for the CENTER reconnection case. All panels share the same 
horizontal axis. Top panel: electron density along the tube for different times in the simulation. The IE- 
positions of the bends (only for the right part of the tube) are shown as circles for each chosen time. Lower 
densities in the central part of the tube correspond to later times in the simulation. The inset figure repeats 
the electron density profile for the time t — 60 s, where each point of the curve corresponds to a grid point 
of the tube, to show that the sub-shocks are well-resolved by the DEFT computer program. Middle panel: 
temperature along the tube for the same times as in the top panel. Large and small crosses indicate the 
beginning (heat front) and end (sub-shock) of the gas-dynamics shocks, respectively. The inset figure repeats 
the temperature profile for the time t — 60 s, where each point of the curve corresponds to a grid point of 
the tube, to show that thermal fronts are well-resolved by the DEFT computer program. The arrow points 
to the location where the back side of the tube (z < 0) becomes tangent to the edge of the current sheet 
(xbottom)- Bottom panel: parallel (along the tube) speed profile for the same times as in the top panel. 
Circles are also the same as in the top panel. Times in seconds are shown for each bend position. 
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Fig. 9. — : Tube arrival at the top of the arcade cartoon. Panel (a): reconnected tube (dark solid line) 
lays on top of the arcade after having descended from the reconnection site. Dashed lines represent the 
underlaying arcade. The legs of the tube connect the parts of the tube that were close to the CS with the 
tube's footpoints at the solar surface. Panel (b): cross section of the tube (solid black circle) before its arrival 
at the top of the arcade (dashed parabolic lines). Panel (c): arcade restoring force. The underlying arcade 
gets compressed and lightly curved due to the arrival of the reconnected tube, and these stronger field lines 
exert an upward force on the tube. 
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Fig. 10. — : DOWNWARD moving tube for the CENTER reconnection case. The two panels share the 
same horizontal axis that corresponds to the right side of the reconnected tube. Top panel: the parallel 
component of the speed for times between 60 and 150 s, at five second intervals. The first time (dashed line) 
corresponds to t = 60 s, which is the last time shown in Figure [7] Time t = 150 s is shown with a thicker 
line. Circles show the x-position of the bends. The "SB" arrows point toward the second bend locations at 
t = 65 s and t = 70 s. The "SB-B" arrow indicates a time when the second bend and the original bend are 
close. Bottom panel: electron density profile for the same times as in the top panel. The arrow points to 
the location where the tube becomes tangent to the edge of the current sheet. 
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Fig. 11. — : DOWNWARD moving tube for the CENTER reconnection case, with the same format as Figure 
IH The time evolution of the mean temperature of the tube (left vertical axis) is shown with a dotted line. 
The total emission measure per flux of the tube (right vertical axis) is also shown for the same times with 
a dashed line. Arrow "C" indicates the time when the center of the tube arrives at the top of the arcade. 
Arrow "B" indicates the time when the bends arrive at the top of the arcade. 
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Fig. 12. — : DOWNWARD moving tube for the TOP reconnection case, , with the same format as Figure 

HU 
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<T> (normalized) 

Fig. 13. — : Normalized emission measure as function of the normalized mean temperature of the tube for 
DOWNWARD moving tubes. The solid line corresponds to the CENTER reconnection case, the dotted line 
to the BOTTOM reconnection case, and the dashed line to the TOP reconnection case. 



